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Quantized Berry phases as local order parameters in t-J models are studied. A texture 
pattern of the local order parameters is topologically stable due to the quantization of non- 
■ Abelian Berry phases defined by low-energy states below a spin gap, which exists in the large 

J/t case with a few electrons. We have confirmed that itinerant singlets in the wide class of 
t-J models carry the nontrivial Berry phase n. In the large J/t case for the one-dimensional 
t-J model, Berry phases are uniformly ir when the number of electrons is N — An + 2, 

V : (n = 0,l,2,...). 
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New quantum numbers have recently been proposed to classify states of matters beyond 
q | the Landau symmetry-breaking description, such as ground state degeneracy, non-Abelian 

Berry's phases and edge excitations. 1-4 Such a nontrivial class of matter is identified as possess- 
ing topological order, which can be useful to characterize quantum liquids; the Haldane spin 
chain, 5 the valence bond solid states, 6 and the spin-Peierls system. 7 Also, the entanglement 
entropy has attracted attention, 8 which provides a new physical insight for the gapped quan- 
turn liquids with a nontrivial topological order thorough the bulk-edge correspondence. 9 ' 10 
Moreover, topological quantities such as the Chern numbers, which have successfully charac- 
terized the quantum Hall states 11 and the generic quantum liquids, 3,4 have the advantage that 
these quantities are quantized, which implies they are topologically protected against small 
perturbations. In addition to the Chern numbers, the Berry phases can be also quantized in 
some situations. When the ground state is gapped and is invariant under some anti-unitary 
operation, the Berry phase, which takes an arbitrary real number in principle, is quantized 
to the two values: a trivial value or a nontrivial value n (mod 2tt). This scheme has been 
proposed and applied for gapped quantum liquids. 12 ' 13 Since the Berry phase is defined by a 
local SU (2) twist at a specific link, the Berry phase can be used as a local order parameter, 
which does not have any classical correspondent since it is gauge dependent. Note that the 
expectation value of any observable is invariant under unitary transformations of bases, which 
means that it is a gauge invariant. In dimerized Heisenberg systems, the Berry phases break 
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the translational symmetry: the Berry phases at strong bonds (or at weak bonds) are ir (or 
zero). It implies that dimer singlets are mostly localized at the strong bonds. The scheme 
can be applied to the strongly correlated electron system as well. 12 To take an advantage of 
topological stability against small perturbations such as randomness, a finite gap above the 
ground state is required. If one does not require the topological stability, a Berry phase can 
be obtained even if there is no finite gap. 14 

In the present work, the scheme to calculate the topologically stable Berry phase for a 
gapped system is applied to the t-J model as one of the strongly correlated systems. 12 ' 13 
In contrast to the localized singlets in the dimerized Heisenberg systems, singlets in the t-J 
model are itinerant. Therefore there can be gapless charge excitations above the ground state, 
which implies existence of a low energy cluster of eigenstates even if the spin excitation is 
gapped. Then, we need to use a non-Abelian Berry connection, which is defined by all the 
eigenstates in the cluster. 12 The Berry phase in this case is defined by taking a trace over 
the eigenstates. It reminds us that the standard order parameter as an expectation value of 
an operator is calculated by taking a trace with the density matrix when the system is in a 
mixed state; besides, when the ground state is degenerate, one needs to take an average over 
the degenerate eigenstates at zero temperature. 

The t-J model is an effective model of the large U repulsive Hubbard model, which de- 
scribes the hole-doped Mott insulator. The study on two dimensional(2D) t-J models was 
stimulated by the discovery of high-T c superconductivity of the copper-oxygen planes. 15 ' 16 
In addition, one-dimensional(lD) t-J models have received much attention recently because 
of the common properties between ID and 2D strongly correlated electron systems 17 and, 
more directly, possibilities of quasi-lD superconductivity. 18 Motivated this proposition, many 
theoretical studies on ID t-J models were focused on superconducting states. For simplicity, 
we focus on the large J/t region in the paper. 

The t-J model is defined on the subspace without double occupancy as 

H = H t + Hj, (1) 
H t = 5>^cJ Cj 7>, (2) 

ij 

H J = ^2 Ji i { Si ' S i J 1 ) ' ^ 

ij 

where cj = (c|^,c||), rij = cjcj, S, = ^ c \ ac i an d <x is a vector form of the Pauli matrices. Here 
the operator V = \\ ri {l — cjt c *t c !| c *l) is the projection to states with no double-occupancy. 
In this paper, we limit ourselves to a 2D square lattice or a ID chain under the periodic 
boundary condition with the lattice size L and the number of up (down) electrons N^(N^). 
Hereafter, the hopping matrix tij has a non-zero value t only on nearest-neighbor links, while 
dimerization or next-nearest-neighbor interactions can be introduced for the spin exchange 
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Jij. We take t = 1 as a unit of energy in the numerical calculations. Let us describe the 
method to calculate Berry phases. For a parameter-dependent Hamiltonian, H = H(9), the 
Berry phase 7 of an M-fold multiplet is customarily given by the U{1) part of non-Abelian 
Berry connection as 27 = j c TrA(0) , where the Berry connection A is the M x M matrix 
defined as {A(6)) mm , = (i/j m (9)\ d\ip m i(9)) = (VVi(0)|"ga|Vw(0)) d# and states are normalized 
eigenstates of the Schrodinger equation; H(9)\ip m (9)) = E m (9)\ip m (9)) , (1 < m < M). 19 
Note that when the M states are degenerate at 9 = a standard order parameter at zero 
temperature is defined as (O) = ^ m (V'm(0)|C|V , m(0))/M and has some analogy with the 
non-Abelian Berry phase. The main differences are the differentiation and the integration for 
the Berry phase, which lead to the gauge dependence. 

The Berry phase 7 is proved to be real and has ambiguity due to gauge freedom of 
eigenstates. To avoid it, we should fix the gauge. Following ref. 3, the gauge-fixed states can 
be obtained from P\(j) m ) with a generic basis set \(j) m ) and the gauge invariant projection 
operator P defined as P = \ip m (9)) {ip m (9)\. However, in this paper, we use another way 
to calculate the Berry phase by introducing a gauge-invariant Berry phase for the lattice 
analogue of the Berry connection. It is defined as follows by discretizing the parameter space 
of 9 into K points: 12 ' 20 

K 

7 = Jim IK, IK = -VargdetC(6» fe ), (4) 

K — >oo L — ' 

fc=l 

where C(9k) is the M x M matrix defined as (C(9k)) mrn / = {ipm(6k)\ipm'(6k+i)) under the 
periodic condition: 9k+i = 9\. The Berry phase in a continuum is gauge dependent but the 
one defined here for the discretized parameter space is gauge invariant. The gauge invar iance 
here means that it is invariant under gauge transformations after fixing the discretization. This 
Berry connection can be considered as the connection of a ID analogue of the lattice gauge 
theory. Note that the periodic condition in the parameter space, \ip m (0K+i)) = IV'm(^i)}) 
guarantees the gauge invariance of the Berry phase here. As a generic parameter 9 in the 
definition of the Berry phase, we use a local spin twist in the present study: 

S+S- + S-S+ e^SfSj+e-MvSrSf, (5) 

where Sf = Sf ± iSf. Since the Berry phase is defined by the local twist at each link ij, the 
Berry phase is used as a local order parameter on the link. Note that 9ij is introduced only 
to the spin exchange terms Hj, and it does not modify the other terms in the Hamiltonian. 
We denote the Hamiltonian with 9{j as H(9ij). 

The gauge-invariant Berry phase here is also quantized when the Hamiltonian H(9ij) is 
invariant under an anti-unitary operator, i.e., [H(9ij),Q] = 0. The anti-unitary operator in 
the present case is the time-reversal operator written as = KU with complex conjugation K 
and a unitary operator defined as U = f\i e^^~ S i \ The quantization of the K-discretized 
Berry phase 7x is proved in the same way as in refs.12 and 13. Since the Berry phase is 
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quantized (7 = or ir) and is used as a link-variable as discussed above, each link has one of 
three labels: "0-bond", "7r-bond", or "undefined (or gapless)". As a whole system, a texture 
pattern of the local order parameters is obtained. 

Before we show the results, we comment on some technical aspects. We can obtain 7 even 
for small K without numerical error due to the quantization of jk- For example, in the one 
singlet case (L = 2, N-\ = N± = 1), K = 3 is enough to obtain the correct value, 7 = tt, in the 

K 

large K limit. On the other hand, absolute value of = det C(6k) is used as a convenient 

k=i 

criterion for the convergence. We have checked it for all the results shown below. 

In the t-J model with a few electrons, the spin gap is finite but the charge excitation is 
gapless. Then, low energy states below the spin gap are treated as an M-fold multiplet to 
calculate the Berry phase. To understand the spin gap in a few electrons case, we consider the 
large J limit. For simplicity, let us consider the one singlet case, = iVj_ = 1 for a while. The 
M-fold degenerate states below the spin gap form a multiplet, where M is equal to the number 
of links and the spin gap is of the order of J. These are spanned by the localized singlet states 
\l), where the link number I indicates the position of the localized singlet (1 < I < M). These 
states are bundled as a multiplet. Then, we consider the hopping process of a singlet as a 
perturbation in a similar way to refs. 16 and 21. In the eigenspace spanned by |/), the effective 
Hamiltonian is obtained by the second perturbation theory as 

H cS = £|0(Z|i^|0<n, H^ = H t¥ ± ir H t . 

iv J 

When we switch on hopping process, the singlet starts to move around with the effective 
hopping t s = y. The degeneracy is lifted by and its bandwidth is of the order of t s . 
When J/t is large enough, the spin gap above the band is stable in the thermodynamic limit, 
L — > 00, as shown in Fig. 1. The charge gap, which is a small gap between the states in the 
multiplet, is of the order oit s /M and becomes zero in the thermodynamic limit. Note that the 
number of the states is M = L for a ID chain or M = 2L 2 for a 2D square lattice. Although 
the degeneracy of the multiplet is macroscopic, the Berry phase of states below the spin gap 
is numerically stable. 

Figure 2(a) shows a texture pattern for one singlet states (iVj = = 1) in the t-J model 
on a 4 x 4 lattice under the periodic boundary condition. When the spin exchange J/t is large, 
a finite gap exists above the M = 32 states. As shown in Fig. 2(a) the uniform 7r-bonds are 
obtained. In addition, when the dimerization of the spin exchange is introduced as J s > J w , 
the dimerization gap opens up. The Berry phase of the states below the dimerization gap will 
be 7r for strong bonds with J s and for weak bonds with J w . For example, as shown in 
Fig. 2(b), when eight links with the strong interaction are alternately distributed along the 
x— axis on a 4 x 4 lattice, the dimerization gap opens above the M = 8 states and the texture 
pattern coincides with the distribution of the strong bonds with J s . 
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Fig. 1. Energy diagram as a function of a local spin twist 9 in the ID t-J model with Nf = N± = 1 
and J/t = 6. The system size is L = 12, 18, 24. 




Fig. 2. A texture pattern of the local order parameters for the one singlet states (JVf = A/j, = 1) 
in the i-J model on a 4 x 4 lattice under the periodic boundary condition. The shadowed links 
denote the 7r-bonds and the others are the 0-bonds. (a) The local order parameters for the M = 32 
states for the t-J model with J/t — 8. (b) The local order parameters for the M = 8 states of the 
dimerized t-J model with the dimerization J s / J w — 2 and J w /t = 8. The strong bonds with J s 
correspond with the 7r-bonds and the other links are weak bonds with J w . 



To explain the results in Fig. 2, let us consider the perturbation from the large J/t limit 
in the ID case for simplicity. In the large J/t limit, the Berry phase of the multiplet is given 
as a sum of respective Berry phases of the localized singlet states \l). Each state gives the 
Berry phase ir only at the specific link where the singlet exists. Then the Berry phases of the 
multiplet are tt uniformly. It does not depend on the dimensionality. 

Unless the spin gap closes, any perturbation cannot modify the texture pattern of the local 
order parameters even if the level-crossing within the multiplet occurs. This is the topological 
stability of the quantized Berry phase. Then, the result at J/t = 8 shown in Fig. 2(a) should 
be the same as uniform n Berry phase in the large J/t limit, because the spin gap does not 
close at J/t = 8. In any dimension, the situation can be the same. In addition, the result 
shown in Fig. 2(b) is the same as in the case J /t = oo and J w /t = 0. 

It is interesting to consider the Berry phase in other t-J models. In addition to the simple 
t-J chain, three modifications can be of interest as spin-gapped systems at finite filling in the 
ID case; (1) the dimerized t-J chain 22 by putting J^i+i = J(l + (— 1) % 6), (2) the t-J- J' chain 23 
by putting J^i+i = J, = J' and (3) the t-J-V chain 24 by adding Hy = V ]T\ mrii+i. We 
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Fig. 3. Energy diagram for the ID t-J model with L = 8, J/t = 8 as a function of J'/t or V/t in (a) 
and S in (b). The solid lines (dotted lines) are energy levels for the local spin twist 6 = (6 = ir). 
When we consider one of three perturbations of J'/t, V/t and <5, other two parameters are fixed 
to zero. The spin gap exists around E = —8 ~ —4. In addition, the dimerization gap opens in (b) 
at the middle of the multiplet. 



have checked numerically that small perturbations of 5, J' , and V do not close the spin gap 
in the iVj = iVj_ = 1 case as shown in Fig. 3. The itinerant singlet gives uniformly tt Berry 
phase as shown in the wide class of the t-J models. It should be noted that the dimerization 
gap opens in the middle of the multiplet in the dimerized t-J chain and the Berry phase of 
states below the dimerization gap shows the simple alternation of 7r's and zeros: 7 = tt (or 
0) on a link with the strong (or weak) exchange. The Berry phase gives clear classification of 
the two topologically different phases for the dimerized t-J chain as well as for the dimerized 
Heisenberg chain. 13 The differences are (i)the Berry phase in the present paper is that for the 
multiplet and (ii)the Berry phases of states below the spin gap are tt uniformly. 

Let us again consider the TVf = > 1 case in the simple t-J chain. In the large J/t limit, 
the system is in the phase-separated phase and all the electrons form a Heisenberg chain 
island. 25 The length of the island is N = N-^ + and the translation of this island generates 
L-fold degenerate states. The Berry phase of each degenerate state is obtained from that in the 
Heisenberg chain. Note that the Berry phases of the Heisenberg chain with even iV sites under 
the open boundary condition turn out to be tt on the left and right boundaries and show the 
simple alternation of 7r's and zeros on the other links as mentioned in the dimerized case. After 
the summation of respective Berry phases of L-fold degenerate states, it follows that the Berry 
phase of the L-fold multiplet in large J/t limit for = N ) = N/2 and L > 2N is uniformly 
tt on every links when N = An + 2. It is uniformly zero when N = 4n, (n = 0, 1, 2, . . .). The 
perturbation with Ht lifts the degeneracy of the large J/t limit. However, the texture pattern 
of the Berry phases is protected until the spin gap closes. This is numerically confirmed for 
various L at fixed N in the low density limit. Note that the spin gap closes at finite filling, 
i.e. infinite iV at fixed N/L, which is the case of usual phase separation in the ID t-J model. 

In conclusion, quantized Berry phase as a local order parameter has been calculated in 
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the t-J models. Comparing the previous study in Heisenberg models, 12 ' 13 there are gapless 
charge excitations below the finite spin gap which exists in the large J/t case. Even if the 
number of states below the gap is large, the Berry phase has been obtained by the exact 
diagonalization method and Eq. (4) successfully. As a result in the t-J models, the Berry 
phase for the itinerant singlet is uniformly ir on nearest-neighbor links in the ID and 2D 
cases. The itinerant singlet carries the Berry phase it in addition to charge, while the singlet 
does not carry spin. Due to the topological stability of the quantized Berry phase, the picture 
of the itinerant singlet obtained in the present study (for large J) has been valid in the wide 
class of t-J models until the spin gap closes. The ID t-J model has the uniform Berry phase 
7r for the L-fold multiplet below the spin gap especially when the number of electrons N is 
4n + 2. 

Moreover, the dimerized t-J model has been classified by a texture pattern of the Berry 
phases in regard to the dimerization gap. In general, the Berry phase can be defined for each 
gap and are protected until the gap closing. This method will be useful even for the frustrated 
electron system. The texture pattern of the Berry phases can be used to find a path to a 
simple strong limit without gap closing as the results shown in the present paper have the 
corresponding strong limits. These strong limits tell us the topological property of the phase 
by using the adiabatic continuity. 
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